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Abstract 



A measurement of the J/tp and tp{2S) production cross sections in pp collisions at 
\/s = 7 TeV with the CMS experiment at the LHC is presented. The data sample cor- 
responds to an integrated luminosity of 37 pb^^. Using a fit to the invariant mass and 
decay length distributions, production cross sections have been measured separately 
for prompt and non-prompt charmonium states, as a function of the meson transverse 
momentum in several rapidity ranges. In addition, cross sections restricted to the ac- 
ceptance of the CMS detector are given, which are not affected by the polarization 
of the charmonium states. The ratio of the differential production cross sections of 
the two states, where systematic uncertainties largely cancel, is also determined. The 
branching fraction of the inclusive B — ?► xp{2S)X decay is extracted from the ratio of 
the non-prompt cross sections to be: 

B{B t/'(2S)X) = (3.08 ± 0.12 (stat.+syst.) ± 0.13 (theor.) ± 0.42 (Spdg)) x lO^^. 
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1 Introduction 

Quarkonium production at hadron colliders provides importar\t tests of calculations in the con- 
text of both perturbative and non-perturbative quantum chromodynamics (QCD), via measure- 
ments of production cross sections and polarizations. 

The ]/tp and ip (2S) mesons can be produced in proton-proton (pp) collisions through two mech- 
anisms: prompt mesons directly produced in the primary interaction and non-prompt mesons 
from the decay of directly produced b hadrons. In addition, ]/ip production can also occur via 
decays of heavier charmonium states, both S-wave (the ^{2S) itself) and P-wave (the three Xc 
states). The determination of the latter contribution is challenging both theoretically and exper- 
imentally, because it requires detection of the low-energy photons from Xc decays. In nonrela- 
tivistic QCD (NRQCD) models HI El, by adding a contribution to prompt charmonium produc- 
tion through colour-octet states ||3l, a satisfactory description of prompt ]/rp and ^^(25) meson 
cross sections at the Tevatron f?] has been obtained. However, in these calculations the frac- 
tion of originating from Xc decays must be assumed from experimental measurements with 
large uncertainties, which makes the tp{2S) mesons cleaner probes of NRQCD predictions. 

Non-prompt ]/tp and ip{2S) production can be directly related to b-hadron production, provid- 
ing a measurement of the b-hadron cross section in pp collisions. Past discrepancies between 
the Tevatron results and the next-to-leading-order (NLO) QCD calculations have been resolved 
using the fixed-order next-to-leading-log (FONLL) approach and updated measurements of the 
h ^]/xp fragmentation and decay IHllll. 

Measurements of the prompt and non-prompt production cross sections of ]/ip mesons decay- 
ing to muon pairs was published using the first Compact Muon Solenoid (CMS) data [7J, as 
well as data from other LHC experiments |l8]- [T0ll . The present analysis extends the CMS result 
with a larger amount of statistically independent data and complements it by providing a mea- 
surement of the ip{2S) production cross section, as well as the ratio of the two cross sections. 
Higher trigger thresholds induced by the increased LHC luminosity do not allow the current 
measurement to reach charmonium transverse-momentum (px) values as low as in Ref. [/J, 
but the high-px reach is increased by the much larger amount of data. The advantage of mea- 
suring the ^{2S) to ]/tp cross-section ratio lies in the cancellation of several experimental and 
theoretical uncertainties. 

The polarizations of the J/i/? and ^{2S) states affect the muon momentum spectrum in the lab- 
oratory frame, thus influencing the charmonium acceptance and, as a consequence, the ex- 
tracted cross sections. Therefore it was decided to present the results in two different ways. 
The first approach assumes unpolarized production for prompt ]/tp and ip{2S) whereas non- 
prompt mesons are assumed to have the polarization generated by the EvtGen Monte Carlo 
program IfTTl , corrected to match the most recent measurements ||T2| . Typical changes of the 
measured cross sections resulting from using hypotheses of full longitudinal or full transverse 
polarizations are also given. The second approach provides results restricted to the phase-space 
region of the CMS muon detector acceptance, in order to avoid corrections which depend on 
the unknown polarizations of the two charmonium states. 

The paper is structured as follows. In Section |2| a brief description of the CMS detector is 
provided. In Section|3| the data and Monte Carlo samples are presented, and the event selection 
is described, while Section [1] presents the method to extract the total ]/rp and t/j(2S) yields. In 
Section |5| corrections for acceptance and efficiency are explained, which are used to determine 
the inclusive cross sections, as discussed in Section |6] Section [7| describes the method to extract 
the J/t/J and ^{2S) non-prompt fractions of the total yields. In Section Isl the results of the 
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differential cross sections, the non-prompt fractions, the cross-section ratios, and the inclusive 
branching fraction i5(B — )■ xp{2S)X) from the ratio of the non-prompt cross sections are given. 

2 The CMS detector 

A detailed description of the detector can be found elsewhere ||T3| . The central feature of the 
CMS apparatus, composed of a central barrel and two endcaps, is a 6 m diameter superconduct- 
ing solenoid producing a 3.8 T magnetic field. Within the magnetic field volume are the silicon 
tracker, the crystal electromagnetic calorimeter, and the brass /scintillator hadron calorimeter. 

The coordinate system adopted by CMS has the origin at the nominal collision point, the y axis 
pointing vertically upward, and the x axis pointing radially toward the centre of the LHC ring. 
The z axis points along the anti-clockwise beam direction defining a right-handed coordinate 
system. The polar angle 9 is measured from the z axis. The pseudorapidity of a particle is 
defined as rj = — ln[tan(0/2)], which approaches the rapidity y = 0.51n[(E + cpz)/ {E — cpz)] 
in the ultra-relativistic limit, where E and pz are the particle's energy and longitudinal momen- 
tum. 

Muons are detected in the pseudorapidity range \rj\ < 2.4 by three types of gas-based detectors 
embedded in the steel return yoke: drift tubes in the barrel, cathode strip chambers in the 
endcaps, and resistive plate chambers in both the barrel and endcaps. 

The silicon tracker consists of the inner pixel-based detector followed by layers of microstrip 
detectors. The strong magnetic field and the good position resolution of the silicon tracker 
enable the transverse momentum of muons matched to reconstructed tracks to be measured 
with a resolution of ~ 1.5 % for pj smaller than 100 GeV/c. 

The first level (LI) of the CMS trigger system, composed of custom hardware processors, uses 
information from the calorimeters and muon detectors to select the most interesting events. 
The high level trigger (HLT) runs on a processor farm to reduce further the rate before data 
storage. 

3 Data selection and event reconstruction 

This analysis is based on a data sample collected in 2010 with the CMS detector, in pp collisions 
at a centre-of-mass energy of 7 TeV. The sample is selected to have consistent trigger require- 
ments for the data used in the analysis and without overlap with the sample used in Ref. [Zl- It 
corresponds to a total integrated luminosity of 36.7 ± 1.5 pb^^ [14]. During this data-taking pe- 
riod, there were on average 2.2 inelastic pp collisions per bunch crossing at the CMS interaction 
region. 

The ]/tp and ^{2S) mesons are reconstructed in the fi^ji^ decay channel. This analysis is based 
on events selected by dimuon triggers that exploit advanced processing at the HLT level. Infor- 
mation from all three muon systems, as well as from the tracker, are used to make the trigger 
decision. Both muons are required to be consistent with a LI muon signal, requiring at least two 
independent segments in the muon chambers, and to be matched to a track reconstructed in a 
region of interest defined by the LI seed. No explicit requirement on the transverse momentum 
Pj is applied. 

Simulated events are used to tune the selection criteria, check the agreement with data, com- 
pute the acceptance, and derive efficiency corrections, as well as for systematic studies. Prompt 
J/t/j and i/'(2S) events are simulated using pythia 6.422 [15J, which generates events based on 



3 



the leading-order colour-singlet and colour-octet mechanisms. Colour-octet states undergo a 
shower evolution. We use the NRQCD matrix element tuning obtained by fitting NRQCD 
calculations to CDF data l|T6l[l7l . In the absence of consistent theoretical and experimental in- 
formation about the j/ip and J/'(2S) polarizations, the dilepton decay distribution is assumed to 
be isotropic. Simulated events with b-hadron decays are also generated with pythi A, and the b 
hadrons are forced to decay inclusively into ]/xp and ^{2S) using the EvtGen package. Photon 
final-state radiation (FSR) is implemented using PHOTOS IIT8l[T9l . 

The off-line event selection, only briefly summarized here, is very similar to the one used 
in Ref. [7J. Muon candidates are reconstructed from the combination of muon-detector and 
silicon-tracker hits. The muons are required to pass the following criteria in the tracker: have 
at least 10 tracker hits, at least two of which are required to be in the pixel layers; have a per 
degree of freedom smaller than 1.8; and pass within a cylinder of radius 3 cm and length 30 cm 
centred at the beam-spot centroid position and parallel to the beam line. 

To select events with ]/xp or xp{2S) decays, muons with opposite charge are paired and their 
invariant mass is computed. The mass is required to be between 2.5 and 4.7GeV/c^. The two 
muon trajectories are refitted with a common vertex constraint, and events are retained if the 
probability of the fit is larger than 1%. If more than one muon pair is found in the event, the 
one with the largest vertex probability is retained. 

The dimuon LI triggers include a veto, whose specific criteria depend on the type of muon 
chamber and the region of the detector: this rejects muon signals whose spatial separation in 
the muon stations is too small, in order to avoid spurious dimuon signatures from a single 
muon. As a consequence, the dimuon sample is split in two, depending on the signed differ- 
ence in azimuthal angle {Acp) between the positively and the negatively charged muons. Muons 
that bend towards each other in the magnetic field are called type-C ("convergent") dimuons, 
while muons that bend away from each other are type-D ("divergent") dimuons. 

Dimuons of type D are much less affected by the trigger veto, while dimuons of type C may 
cross at the muon stations. This causes sizeable correlations between the two muon detection 
efficiencies, and this effect is larger in the forward region. Therefore, in addition to the above 
requirements, all type-C dimuons are rejected for the inclusive cross-section measurements. 
This corresponds to a 48% reduction in the yield. For the non-prompt fraction determination 
(which is largely efficiency independent), type-C dimuons are only rejected in the high dimuon 
rapidity region 1.6 < \y\ < 2.4. 

The momentum measurement of charged tracks in the CMS detector has systematic uncer- 
tainties that are due to imperfect knowledge of the magnetic field, modelling of the detector 
material, and sub-detector misalignment. These effects can shift and change the width of the 
mass peaks of dimuon resonances. In addition to calibrations already applied to the data 1201 - 
[22], residual effects are determined by studying the dependence of the reconstructed dimuon 
peak shapes on the muon kinematics, as was done in Ref. ||Z|- 

Because of the large difference in the branching fractions to dimuons of the two states, the 
measured ^{2S) yield is much smaller than the ]/xp yield. For this reason, different binnings 
are used for the differential cross sections in pj and \y\. 

4 Inclusive yield determination 

Two methods are used to extract the inclusive yields from the fi^ji^ invariant mass distribution, 
either fitting the ]/^p peak alone in a restricted mass window, or fitting the combined J/i/' and 
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tp{2S) distribution. Yields are derived using an extended unbirmed maximum-likelihood meth- 
od. 

In both types of fits, the sum of a Gaussian and a Crystal Ball p3| function is used for the de- 
scription of the signal, simultaneously taking into account FSR and rapidity-dependent resolu- 
tion variations. In the J/t/^-only fits, an exponential function is used to describe the background. 
Figure [l] (top) shows an example of a fitted mass distribution. The mass resolution ranges from 
about 20 MeV/ in the low rapidity region to 35 MeV/ at intermediate rapidities, up to about 
50 MeV/ in the forward, high rapidity region. 

In the second type of fits, the two mass peaks and the background are fitted simultaneously. 
For the J/t/; and 1/^(25) signal peaks, the same probability density function (pdf) is used, with 
the following constraints on the parameters: the ratio of the central values of the two masses 
is fixed to the world average value ||24| : the widths, scaled by the nominal mass values, are 
constrained to be the same, and the parameters describing the asymmetric tail of the Crystal 
Ball function are constrained to be equal. The background is modelled by two exponentials. 
Figure [T] (bottom) shows an example of a fitted mass distribution. 
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Figure 1: Top: The ji^ji invariant mass distribution in the ]/xp region and the result of the fit 
for the bin: 0.9 < \y\ < 1.2, 12 < pj < 15 GeV/c. Bottom: ]/xp and xp{2S) mass distribution fit 
for the bin: < 1.2, 8 < pj < 9GeV/c. The solid and dashed lines represent the total fits and 
their background components, respectively. 
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The systematic uncertainties on the mass distribution fits are estimated by changing the ana- 
lytical form of the signal and background pdf hypotheses (a single Crystal Ball function is used 
for the signal and polynomial pdfs for the background) in the two types of fits. The largest 
variation in the yield of each fit is taken as the systematic uncertainty. 



5 Acceptance and efficiency 

As discussed in Section [l| measurements will be presented using two different approaches. In 
the first approach, the observed number of]/xp events is corrected for the detector acceptance 
and reconstruction efficiency in every bin in which the cross section is measured. As the ac- 
ceptance is strongly dependent on the assumed polarization of the charmonium state, in the 
second approach we provide measurements exclusively within the CMS detector acceptance, 
where only detector efficiency corrections are made and without any polarization-related un- 
certainties. 



5.1 Acceptance 

The acceptance reflects the geometrical coverage of the CMS detector and the kinematic reach 
of the muon trigger and reconstruction, constrained by the amount of material in front of the 
muon detectors and by the track curvature in the magnetic field. 

In the simulation, both muons are required to be within the geometric acceptance of the muon 
detectors. A single muon is defined as detectable if it satisfies the following requirements at 
generator level: 

p^>4.0GeV/c for \rj^\<l.l 

> 3.3 GeV/c for 1.2 < |?/^| < 2.4 . (1) 



The ]/-tp acceptance A is defined as the fraction of detectable ]/\p —> ji^ji decays, as a function 
of the generated dimuon transverse momentum pj and rapidity y, 



A[pT,y,A0} - — , , . , (2) 

Ngen(PT/ y, Ae) 

where N^^^ is the number of detectable ]/ip events in a given {pj, y) bin, and Ngen is the cor- 
responding total number of generated ]/xp events in the Monte Carlo (MC) simulation. An 
analogous definition holds for xp(2S). The parameter Aq reflects the fact that the acceptance 
is computed for various polarization scenarios, which lead to different muon spectra in the 
laboratory frame. 

For the acceptance calculation, a dedicated sample of generated events is used, with no restric- 
tions on the phase space. The large number of simulated events allows a much smaller bin size 
for determining A with respect to that used for the cross-section determination in data. 

To study the effect of ]/xp and tp{2S) polarization on the acceptance, these events are reweighted 
depending on the values of the polar and azimuthal angles as computed in two different frames 
(helicity and Collins-Soper 1125 1 ). The angular distribution for the decay of a / = 1 state into 
fermions is used, which is a function of three independent parameters Ag, A^, and Ag^p: 

3 

W(cos0,(*) = — ; • (1 + Aecos^0 + Ad,sin^0cos2(* + A0rf,sin20cosrf)) . (3) 

2[3 + Ag) 
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The choice of zero for all A parameters corresponds to an unpolarized decay, while Ag = —1 
and \q = +1 correspond to fully longitudinal and fully transverse polarizations, respectively. 

By default, the prompt ]/rp and ip{2S) are assumed to be unpolarized, while the non-prompt 
mesons are assumed to be polarized as generated by EvtGen and corrected to match recent 
measurements, as mentioned in the Introduction. Typical changes of the measured cross sec- 
tions when using alternative polarization scenarios are provided in Section |8j 

Several sources of systematic uncertainty on the acceptance have been investigated: 

• Kinematic spectra. Different px and y spectra of the generated ]/xp and ip{2S) might 
produce different acceptances, as the acceptance is defined by single-muon crite- 
ria. Spectra from theoretical predictions presented in Section |8] have been used to 
recompute the acceptance, and the difference from that obtained with the pythia 
spectrum has been taken as a systematic uncertainty. 

• Final-state radiation. The generated dimuon momentum may differ from the j/ip and 
i/'(2S) momentum, because of FSR. The difference between the acceptance computed 
using the dimuon or the charmonium variables in Eq. (|2j| is taken as a systematic 
uncertainty. 

• B polarization. The ]/tp and ^{2S) mesons produced in b-hadron decays have a dif- 
ferent acceptance with respect to the prompt ones. The corresponding systematic 
uncertainty is evaluated by taking the difference of the default choice (corrected to 
match the experimental results of Ref. [12]) with respect to the one predicted by 
EvtGen. 

• Pj calibration. The muon transverse momenta in data have been calibrated as de- 
scribed in Ref. [7J. A difference in the momentum resolution between data and 
simulated events would also give a systematic uncertainty on the acceptance. The 
acceptance has been computed with simulated muon momenta smeared according 
to the resolution measured with data |26J. The uncertainty on the measured res- 
olution was used to apply an additional smearing on the simulated momenta; the 
acceptance has been recalculated and the shift taken as a systematic uncertainty. 

5.2 Muon efficiency 

The single-muon efficiency is measured from data for muons in the acceptance, as described in 
Refs. Il27l |28| , and is based on the tag-and-probe (T&P) method. For this purpose, independent 
sets of triggers are used for which online requirements either on the muon or the tracker tracks 
are not applied, thus yielding samples which are unbiased with respect to the corresponding 
selections. 

The combined trigger and offline reconstruction efficiency for a single muon is defined as: 



^(m) - ^trig I off ■ ^off I ID ■ ^ID I track ' ^track' (4) 

where ^track offline tracking efficiency, e^Q | track refers to the muon identification in the 

muon systems for a tracker-reconstructed muon, e^^^ | refers to the specific quality require- 
ments applied to reconstructed muons, and e^j.jg | is the probability for an offline recon- 
structed muon to have also fired the trigger. 

The muon identification and trigger efficiencies (e^-j-jg i gff/ ^off i id and ejp i ^rack^ have the 
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strongest Pj and |?/''| dependence and are determined in 15 bins of pj (3.3 < pj < 50GeV/c) 
and 14 bins of Irjf'l (0 < < 2.4), allowing an adequate description of the turn-on efficiency 
curves. Since the tracking efficiency is almost constant for this momentum and rapidity range, 
broader bins are used. 

The efficiency to detect a dimuon event is expressed as: 

= e{^^) ■ e{}i-) ■ p ■ evertex. (5) 

where £{}i^) and are the single-muon efficiencies, and %ertex is the efficiency of the 

vertex requirement, calculated from the data by determining the yields in regions of large 
Pj and \y\hy alternatively applying and not applying this requirement. The p factor, defined 
by Eq. (jSjl, represents a correction to the efficiency factorization hypothesis: it accounts for the 
finite size of the (pj,?/^) bins and, more importantly, for the possible bias introduced by the 
T&P measurement, due to correlation effects as discussed in Section |3] In order to determine p, 
the efficiencies have also been evaluated using T&P techniques on simulated events and their 
product has been compared with the true dimuon efficiency. Except for some bins at high pj, 
the values are found to satisfy |1 — |0| < 10%. 

For the acceptance-corrected cross-section results, the acceptance and efficiencies are combined 
into a single factor, which is computed for each {pj, y) bin and is defined as: 



where the average is taken over the data events in each bin, using the "fine-grained" bins of the 
acceptance and the event-by-event efficiency obtained from the single-muon efficiencies using 
Eq. (|). 

Similarly, for the results which are not corrected for acceptance, the efficiency factor is deter- 
mined as: 



Ne, 



Two sources of systematic uncertainty in the efficiency are considered: 

• The uncertainties on the measured muon efficiencies propagate as systematic errors 
on the cross section measurement through the correction factor (^)bin (or (^)bm)- 
The effect has been estimated on a statistical basis in each bin by performing Monte 
Carlo pseudo-experiments, in which the muon efficiencies were varied randomly 
according to a probability density built by joining the left and the right side of two 
Gaussians with different widths, in order to allow for asymmetric errors. The r.m.s. 
of these correction factors in each bin has been taken as the systematic uncertainty 
associated with the single-muon efficiency. 

• The full difference 1 1 — p | in each bin is taken as a systematic uncertainty due to the 
efficiency correlation. 
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The inclusive double-differential cross section is given by: 



d2 



{]/lp)■B{]/^p^^+^-) 



Tvfcorr 



J Ldt ■ Apx ■ Ay 



(8) 



where J L dHs the integrated luminosity, Apj and Ay are the pj and y bin widths, S(J/t/^ — ?► 
fi^fi^) is the decay branching fraction of the ]/tp into two muons, and Nj/^"(pT, \y\) is the cor- 
rected yield in a given (pj, \y\) bin. The corrected yield is obtained from the fitted signal 
yield Nj/^ via Nj°^^ = Nj/^ ■ {^)bin in the case where the ]/ip yields are corrected for accep- 
tance and efficiency, and Nj°^ = Nj/^ • (^)bm in the case where the results are uncorrected for 
acceptance. An analogous formula applies for the tp{2S) double-differential cross section. 

Figure |2] shows the measured (fully corrected) J/xp inclusive cross section as a function of pj for 
the various rapidity bins. They are compared with our previous results published in Ref. [TJ, 
which are statistically independent and remain of interest since they partially overlap with the 
present results and cover a lower pj range. A good agreement is observed. In this figure, as 
well as in the cross-section plots of Section |8| multiplicative factors - appearing as additive 
offsets on the log scale - are used to achieve a convenient graphical separation of the measure- 
ments from different rapidity bins. 
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Figure 2: Measured differential cross section for ]/xp inclusive production as a function of pj for 
five rapidity bins, fully corrected for acceptance and efficiency. Also plotted are the results pub- 
lished in Ref. [7J, which extend to a lower pj range. The error bars on data points include all the 
statistical and systematic contributions except luminosity and polarization. The measurements 
have been offset by the numerical values given in the legend for easier viewing. 
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7 Prompt and non-prompt fractions 

To estimate the ]/ip fraction from b-hadron decays, a two-dimer\sior\al fit is performed, in which 
the pdfs and fit procedure are the same as those described in Ref. HZl- The variables used for 
the two-dimensional fits are the dimuon invariant mass and the "pseudo proper decay length" 
£y^, defined as the most probable value of the transverse distance between the dimuon vertex 
and the primary vertex, corrected by the transverse Lorentz boost of the ]/tp. As in Ref. [7j, the 
primary vertex is chosen as the one closest to the dimuon vertex in the z direction. 

The resolution of the pseudo proper decay length is described by a function depending on 
an event-by-event uncertainty determined from the covariance matrices of the primary and 
secondary vertex fits. The uncertainty is used as the r.m.s. of the resolution Gaussian function 
that describes the core of the resolution, while a second Gaussian function with a small relative 
normalization (usually < 1%) parametrizes the effect of incorrect primary vertex assignments. 

The pdf F{ij/^, m^^, a^) for the ]/ip is then: 

F(£j/^, m^,^, ac) = /sig ■ Dsi^{cri>) ■ Fsig(£j/^, ai) ■ Msig(m,,,,) + 

(1 - /sig) • ^Bkg{(^i) ■ fBkg(%^, (^t) ■ MBkg(m^;,), (9) 

where: 

myp(^d = L^n^v^) ^ - iy^\F,sm). (10) 

and k = {Sig, Bkg}. In the equations above: 

• ^Sig(^f/f/) arid MBkg(w2fi^) are the mass pdfs determined for the signal and back- 
ground in Section |4| and /sig is the fraction of signal events in the entire range of the 
fit; 

• fgig'^(%i/)) and pQ^gi^j/tp) are the functional forms describing the £y^ distribution for 
the signal and background, respectively. The signal part is given by the sum of 
prompt and non-prompt components: fsig''(%i/') = /b • fb{^]/^) + (1 - /b) ■ Fpi^y^p), 
where /b is the fraction of }/ip from b-hadron decays, and fp(£j/^) and fb{^]/Tp) are 
the distributions for prompt and non-prompt J/i/^, respectively. The £y^ pdfs for 
prompt signal and background are the same as in Ref. IZl. The non-prompt life- 
time function is described by an exponential decay of the b hadron, with a Gaussian 
smearing function that accounts for the difference between the measured pseudo 
proper decay length and the proper decay length of the b hadron; 

• (7i is the per-event uncertainty of the decay length and Dsig{crf ) and DBkg(c^^) are its 
distributions separately for signal and background [29] . They are obtained from the 
signal region of the invariant mass distribution, after a sideband subtraction, and 
the sideband regions, respectively; 

• Ri and R2 represent the core and tail decay-length resolution Gaussian functions: fi 
is their common mean and s, represent scale factors for the per-event uncertainty, 
which are both left free in the fit to account for initial assumptions on the uncer- 
tainties of track parameters. These functions are convolved with Fl^^'^{£y^) to obtain 
the observed Fk{£]/ip) distributions, including the experimental resolution {k = {Sig, 
Bkg}). 

The background is fitted using the events in mass sidebands and the result is used to fix lifetime 
parameters of the overall fit in the entire mass region. The mass sideband region is defined as 
[2.50, 2.85] and [3.25, 3.35] GeV/c^. 
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For the determination of the ip{2S) non-prompt fraction, the quantity ^^(2S)/ defined as for the 
]/ip case, is computed. In order to constrain the fit and avoid problems due to Hmited statistical 
accuracy, the ]/tp and 1/^(25) samples are fitted simultaneously using the same binning as for the 
i/?(2S) cross-section determination. The lifetime resolution functions Ri and Ri are constrained 
to be described by the same parameters (mean value and scale factors) and the backgrounds to 
have the same fractions of long-lived components. 

The invariant mass sideband regions used for the determination of the background parameters 
are defined as above for the ]/xp, and as [3.35, 3.45] and [3.85, 4.20] GeV/c^ for the t/;(2S). 

Figure [3| shows two examples of the and ^^(2S) distributions with projections of two-dimen- 
sional nts on these dimensions, as well as the prompt and non-prompt components obtained 
as described above. The lower plots in Figure |3] give the pull distributions from the fits, and 
show no systematic structure. 




Ijiv ("ini) V(2S) ('Tim) 



Figure 3: Left: Projection of a ]/tp two-dimensional fit on the iy^ dimension in the bin \y\ < 0.9, 
8 < Pt < 9GeV/c and in the whole mass region [2.50, 3.35] GeV/c^. Right: Projection of a 
]/tp-xp{2S) two-dimensional fit on the ^,p(2S) dimension in the bin: 1.6 < \y\ < 2.4, 12 < px < 
15GeV/c, in the tp{2S) mass region [3.35, 4.20] GeV/c^. The solid lines represent the total fits; 
the prompt, non-prompt and background components are also shown using green dash-dotted, 
red dashed and blue dotted lines, respectively. The fit pull plots show no systematic structures. 

Several sources of systematic uncertainty have been addressed, using mostly the same proce- 
dures as in Ref . [7] . The main additional systematic effect comes from attempting to choose the 
correct primary vertex of the interaction in the presence of pile-up. The sources of systematic 
uncertainty include the following: 

• Primary vertex assignment. In order to estimate the possible effect of pile-up on the 
primary vertex estimation, the primary vertex associated to the dimuon is chosen as 
the one with the largest track ^ pj^, instead of the one closest in z to the dimuon ver- 
tex. The difference between the fitted non-prompt fractions in these two approaches 
is taken as the systematic uncertainty. 
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• Residual misalignment in the tracker. The effect of uncertainties in the measured mis- 
alignment of the tracker modules is estimated by reconstructing the data using dif- 
ferent sets of alignment constants. The largest difference in the fit results with respect 
to the nominal case is taken as the systematic uncertainty. 

• h-hadron lifetime model. An alternative fit method is used, namely the b-hadron life- 
time model used in Ref. [7\, which is based on MC templates; the difference in the 
fitted non-prompt fraction is taken as the systematic uncertainty. 

• Background fits. The effect of a ±100 MeV/c^ variation in the lower limit of the low- 
mass side (upper limit of the high-mass side) of the ]/xp {tp{2S)) sideband boundaries 
is taken as the systematic uncertainty. 

• Resolution model. The nominal (double-Gaussian) model for the pseudo proper de- 
cay length per-event resolution is compared with a model using a single-Gaussian 
shape. The difference is taken as the systematic uncertainty. 

• Different prompt and non-prompt efficiencies. The MC simulation predicts slight dif- 
ferences between the prompt and non-prompt J/t/? and tp{2S) efficiencies, mostly 
because of the different track densities from fragmentation products around the 
muons. These are taken into account; the relative difference is propagated to the 
non-prompt fraction, and taken as the systematic uncertainty. 

Non-prompt fraction results are given in Section |8] and a summary of all the systematic uncer- 
tainties is given in Tables [T] and |2] for the ]/ip and 1/^(25), respectively. 

8 Results 

The prompt and non-prompt double differential cross sections for the two charmonium states 
are obtained by multiplying the measured inclusive cross sections with the fraction of prompt 
and non-prompt events, respectively. In addition the cross-section ratio of the two charmonium 
states is calculated. 

Statistical uncertainties and contributions from the investigated sources to the total systematic 
uncertainties on these cross sections are summarized in Tables[l]and|2] The largest uncertainties 
are due to the efficiency correlations; FSR estimation has a sizeable effect only in bins close to 
the edges of the acceptance. 

8.1 Prompt and non-prompt cross sections corrected for acceptance 

Figures |4] and |5] show the measured prompt and non-prompt cross sections for the ^/xp and the 
i/?(2S) as a function of pj, for the various rapidity bins and corrected for detector acceptance. 
They are compared with theoretical predictions from NRQCD |3| and from FONLL ||5H6l for 
the prompt and non-prompt cases, respectively. Numerical values are also reported in Ref. II30I . 

The NRQCD prediction includes non-prompt production in the ]/xp case caused by feed-down 
decays from heavier charmonia, and can therefore be directly compared with the data. Good 
agreement is found in both the ]/xp and the ^(2S) cases. For non-prompt production the mea- 
surements lie systematically below the FONLL predictions, possibly because of the large un- 
certainty on the B — )■ xp{2S)X branching fraction. In general, for both states, the observed 
differential cross sections seem to fall more rapidly than the FONLL prediction at high pj, and 
this effect is more evident for the J/i/? because of the higher pj reach. 

The NRQCD theoretical uncertainties include those on the feed-down contributions and on 
the colour-octet, long-distance matrix elements determined from fits to the Tevatron data. The 



12 



8 Results 



Table 1: Summary of tfie relative statistical and systematic uncertainties on the non-prompt 
]/tp cross section (in %). The variation over the different px bins is given for the five rapid- 
ity regions. Uncertainties on the prompt cross section are identical, with the exception of the 
non-prompt fraction, where they must be regarded as relative to (1 — /b) rather than to f^. 
Acceptance uncertainties on the FSR are given, excluding the lowest-px bin in every rapidity 
region, where it can be as large as 19% because of acceptance edge effects. 

Ij/I range 0-0.9 0.9 - 1.2 1.2-1.6 1.6-2.1 2.1-2.4 

Quantity Source Relative uncertainty (in %) 

affected 



All cross sections 



m^i^ fits 


Statistical 


1.2 


-8.9 


1.5 


-7.1 


1.6- 


8.4 


1.2 


-3.2 


2.3 


-3.9 


^j/^ fits 


Statistical 


1.0 


-5.9 


1.4 


-4.7 


1.4- 


7.6 


2.1 


-8.3 


4.4 


-7.1 


Efficiency 


Single-muon efficiency 


0.3 


-0.9 


0.2 


-1.6 


0.1 - 


1.4 


0.2 


-1.0 


0.6 


-1.4 




p factor 


1.9- 


-23.2 


1.2 


-7.6 


0.7- 


5.7 


0.8 


-5.4 


3.7 


-6.8 


Yields 


Fit functions 


0.6 


-3.4 


0.4 


-2.8 


0.5- 


2.8 


0.8 


-2.2 


1.0 


-4.2 


Luminosity 


Luminosity 


4.0 


4.0 


4.0 


4.0 


4.0 


Non-prompt 


Tracker misalignment 


0.1 


-2.1 


0.1 


-0.8 


0.0- 


1.5 


0.2 


-3.2 


0.2 


-5.1 


fraction 


b-lifetime model 


0.1 


-3.0 


0.1 


-3.4 


0.1 - 


3.7 


0.2 


-2.6 


0.2 


-6.6 




Vertex estimation 


0.1 


-0.7 


0.7 


-3.0 


0.4- 


3.7 


1.5 


-4.6 


2.3 


-5.0 




Background fit 


0.0 


-0.2 


0.1 


- 1.4 


0.1 - 


1.0 


0.0 


-2.5 


0.1 


- 1.2 




Resolution model 


0.2 


-3.5 


0.0 


-4.2 


0.8- 


3.5 


1.1 


-5.0 


1.1 


-4.4 




Efficiency 


0.4 


-2.1 


0.9 


-3.3 


0.5- 


9.9 


0.3 


-3.3 


1.6- 


-10.5 




Only acceptance- 


corrected cross sections 












Acceptance 


FSR 


0.0 


-1.5 


0.0 


-2.5 


0.0- 


4.2 


0.7 


-8.0 


0.5 


-3.5 


px calibration 


0.0 


-0.6 


0.0 


-0.6 


0.0- 


0.8 


0.1 


-0.6 


0.0 


-0.8 




Kinematic spectra 


0.0 


-0.3 


0.0 


-0.7 


0.0- 


0.7 


0.7 


-3.8 


0.4 


-5.3 




B polarization 


0.0 


-0.5 


0.0 


-0.4 


0.0- 


0.5 


0.1 


-0.8 


0.3 


-1.3 



FONLL theoretical errors include uncertainties on B{B ]/xp X) and B{B — > ip{2S)X), renor- 
malization and factorization scales, b-quark and c-quark masses, parton distribution functions, 
and bfragmentation parameters. 

However, uncertainties on the B — ?► charmonium decay spectrum were not included in the orig- 
inal FONLL prediction. To estimate those, we make use of the EvtGen MC generator, which 
describes B — ?► charmonium decays using a sum of many exclusive modes. We split the decay 
modes into two categories, "high-Q" and "low-Q", if the value of Q in the decay is respectively 
greater than or less than 1.2GeV/c^, where Q = — J^jtrii and the index i runs over the B 
decay products. As low-Q (high-Q) modes yield charmonia with smaller (larger) momentum 
in the Brest frame, they populate different regions of the B-decay spectrum. Two sets of non- 
prompt charmonium MC events are generated according to the following criteria. In the first, 
each high-Q mode branching fraction is increased by its world-average uncertainty [24J or by 
100% of its value if the branching fraction is not measured. Low-Q mode branching fractions 
are decreased by a similar amount, rescaling the sum to unity after this procedure. In the sec- 
ond, the treatment of the high- and low-Q modes is interchanged. The maximum difference in 
the resulting spectra in the two cases is added to the theoretical FONLL uncertainty. 

To investigate the effect of the assumed ]/xp and tp{2S) polarizations on the prompt cross section, 
the acceptance is recomputed for four extreme polarization scenarios corresponding to fully 
longitudinal or fully transverse polarization in the helicity and Collins-Soper frames |25|. This 
produces relative cross-section shifts across the entire kinematic range of up to 18-20% (20- 
25%) for the ]/ip {tp{2S)) in the helicity frame, and 6-15% for both states in the Collins-Soper 
frame. Detailed results can be found in Ref. ||30| . 



8.2 Prompt and non-prompt cross sections uncorrected for acceptance 
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Table 2: Summary of the relative statistical and systematic uncertainties in the non-prompt 
i/?(2S) cross section (in %). The variation over the different px bins is given for the three ra- 
pidity regions. Uncertainties on the prompt cross section are identical, with the exception of 
the non-prompt fraction, where they must be regarded as relative to (1 — /b) rather than to f^. 
Acceptance uncertainties on the FSR are given excluding the lowest- px bin in every rapidity 
region, where it can be as large as 29% because of acceptance edge effects. 

\y\ range 0-1.2 1.2-1.6 1.6-2.4 

Quantity Source Relative uncertainty (in %) 

affected 



All cross sections 



m^ifi fits 


Statistical 


5.6- 


- 14.8 


7.5- 


31.7 


7.3- 


-24.1 


^il>(2S) fits 


Statistical 


4.3- 


-12.7 


5.9- 


38.0 


9.1- 


-26.4 


Efficiency 


Single-muon efficiency 


0.1 


-0.5 


0.1 - 


-0.6 


0.2 


-0.9 




jO factor 


0.7- 


- 13.1 


2.1 - 


-6.6 


2.3 


-9.8 


Yields 


Fit functions 


1.2 


-3.7 


0.6- 


12.1 


3.1- 


-10.0 


Luminosity 


Luminosity 


4.0 


4.0 


4.0 


Non-prompt 


Tracker misalignment 


0.3 


-2.6 


1.5- 


-7.1 


1.8- 


-11.1 


fraction 


b-lifetime model 


0.0 


-2.5 


0.4- 


-7.6 


0.0 


-2.9 




Vertex estimation 


0.0 


- 1.7 


0.2- 


-3.5 


1.2 


-4.2 




Background fit 


1.0 


-6.8 


2.2 - 


10.0 


2.5- 


-15.3 




Resolution model 


0.5 


-3.5 


0.1 - 


-4.6 


0.9- 


-24.9 




Efficiency 


0.5 


-7.8 


0.9- 


-6.3 


0.5- 


-13.8 


Only acceptance-corrected cross sections 


Acceptance 


FSR 


0.0 


-3.9 


0.5- 


-3.4 


0.3 


-4.1 


Pj calibration 


0.2 


-0.5 


0.3- 


-0.5 


0.3 


-0.5 




Kinematic spectra 


0.1 


- 1.2 


0.0- 


-0.9 


0.7 


-2.0 




B polarization 


0.1 


-0.8 


0.0- 


-0.6 


0.2 


- 1.7 



8.2 Prompt and non-prompt cross sections uncorrected for acceptance 

As discussed previously, since the polarization effects are large compared to the measurement 
uncertainties, cross-section values are also reported that are restricted to the CMS muon accep- 
tance region, to allow future measurements of the ]/tp and i/'(2S) polarization to be exploited. 

Figures |6] and |7| show the measured prompt and non-prompt cross sections for the j/xp and the 
i/7(2S) as a function of pj for the various rapidity bins and uncorrected for detector acceptance. 
Numerical values can be found in Ref. II30I . 

8.3 Non-prompt fractions 

The measured non-prompt fractions for ]/xp and ^{2S) mesons, extracted as described in Sec- 
tion[7|and uncorrected for acceptance, are reported in Ref. |30| and shown in Fig.|8] The uncer- 
tainties shown are statistical and systematic, and the measured values are plotted as a function 
of Pj in three rapidity ranges. In agreement with previous measurements HI [7|, we observe 
similar sizes of non-prompt fractions for ]/xp and xp{2S), and an increasing trend with pj. Ac- 
ceptance corrections do not induce significant changes in the non-prompt fractions within their 
uncertainties. 

8.4 Cross-section ratio 

Most of the systematic uncertainties on the acceptances and efficiencies listed in Tables [T] and |2| 
as well as the luminosity uncertainty, cancel partially or fully in the ratio of the i/'(2S) to ]/xp 
cross sections. For this reason we also present the ratio of the two differential cross sections: 
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Figure 4: Measured differential cross section for prompt ]/xp and ip{2S) production (left and 
right, respectively) as a function of pj for different rapidity bins. The error bars on the data 
points include all the statistical and systematic contributions except luminosity and polariza- 
tion. The measurements have been offset by the numerical values given in the legend for easier 
viewing. The coloured (dark) bands indicate the theoretical predictions from NRQCD calcula- 
tions. The lines are added only for illustrative purposes. 
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Figure 5: Measured differential cross section for non-prompt ]/ip and t/'(2S) production (left 
and right, respectively) as a function of pj for different rapidity bins. The error bars on the data 
points include all the statistical and systematic contributions except luminosity. The measure- 
ments have been offset by the numerical values given in the legend for easier viewing. The 
coloured (dark) bands indicate the theoretical predictions from FONLL calculations. The lines 
are added only for illustrative purposes. 
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Figure 6: Measured differential cross section for prompt ]/ip (left) and t/'(2S) (right) production 
as a function of pj for the different rapidity bins. The error bars on data points include all the 
statistical and systematic contributions except luminosity. The measurements have been offset 
by the numerical values given in the legend for easier viewing. The results are not corrected 
for the muon acceptance. 
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Figure 7: Measured differential cross section for non-prompt ]/xp (left) and i/'(2S) (right) pro- 
duction as a function of pj for the different rapidity bins. The error bars on data points include 
all the statistical and systematic contributions except luminosity. The measurements have been 
offset by the numerical values given in the legend for easier viewing. The result is not corrected 
for the muon acceptance. 
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Figure 8: Fitted ]/tp and tp{2S) non-prompt fractions plotted as a function of pj for three rapidity 
regions: < \y\ < 1.2 (top); 1.2 < \y\ < 1.6 (middle); 1.6 < \y\ < 2.4 (bottom). The inner 
error bars represent the statistical uncertainties only, the outer ones are the quadratic sum of 
statistical and systematic uncertainties. 



I 4^W2S))-g(V.(2S)^;.>-) _ N;°--s)(pT.|y|) 

where the ratio R is computed in bins of px and rapidity and the binning is the same as used 
for the tp{2S) cross section. 

The statistical uncertainties affecting R are extracted directly from the simultaneous invariant 
mass fits. The systematic uncertainties are estimated by considering the same sources as for the 
cross sections (except the luminosity and single-muon efficiency which cancel out in the ratio) 
and evaluating directly the variation of the ratio, in order to take correlations into account. 

No significant dependence of R on rapidity is observed; the ratios over the entire rapidity range 
are therefore computed. The resulting prompt and non-prompt cross-section ratios are shown 



8.5 Inclusive B xp{2S)X branching fraction 
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Figure 9: Measured value of R, the 1/^(25) to ]/xp differential cross-section ratio defined in 



Eq. 11 for prompt (left) and non-prompt (right) production, averaged over rapidity and plot- 
ted as a function of pj. The left plot also includes the comparison with the NRQCD predic- 
tion, while the right plot shows the predictions of the theoretical models used to determine 
B{B — > tp{2S)X), after the latter have been rescaled to the fitted value given in Eq. (12 1. The 
shaded bands show the uncertainties on the theoretical predictions. The error bars give the 
total uncertainties on the measurements; polarization uncertainties are not included. 

in Fig.|9]as a function of pj. Numerical values of rapidity-dependent and integrated ratios are 
given in Ref. ||30| . 

The assumptions on polarization also affect the prompt cross-section ratio measurement. In 
a plausible scenario [31], the polarizations of the directly produced ]/xp and t/j(2S) states are 
assumed to be the same. Therefore the uncertainty on the ratio comes only from the difference 
between the polarization of the directly produced mesons and the polarization of the ]/xp com- 
ing from decays of P-wave states (Xci and Xci), for which the maximum possible variations are 
considered. Using the measured feed-down fractions measured at CDF ||32j [33, this leads to 
the definition of the two extreme scenarios: 

. Af ^ 1^ ^ 0.445; 



• A 



-1, A^^^ = -0.647; 



which result in changes to the measured prompt cross-section ratio by 12 — 20%. 

8.5 Inclusive B tp {2S)X branching fraction 

The non-prompt ip{2S) cross-section results can be used to determine B{B tp{2S)X), the 
average inclusive branching fraction of all weakly decaying particles containing a b quark to 
ip{2S). 

Since the results are determined only for a limited range of phase space, a theoretical assump- 
tion is needed to extrapolate to the full phase space. The most precise result is obtained using 
the non-prompt cross-section ratio, where most theoretical uncertainties cancel. The FONLL 
model is used for the result, taking as an alternative the EvtGen prediction to determine a 
systematic uncertainty. 

For both models, the predicted ratio is computed for each pj bin used in the measurement. 
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9 Summary 



assuming the world-average values, listed by the Particle Data Group (PDG) 1241 , of ^Bpdg ~^ 
]/xp X), BpijG{xp{2S) —> }i^}L^), and SpdgQ/'/' ~^ H^H ) (the branching fractions taken from 
the PDG are indicated as Spdg)- In both models the production cross section for each type 
of B meson or baryon is weighted using the values of the fractions /(b — > B"), /(b — > B^), 
/(b —7- Bg ), and /(b — > A^), taken from LEP and Tevatron measurements. The predictions are 
then fitted to the data points, leaving only the normalization (N) as a free parameter. A good 
agreement in the shape of the pj distribution is found for both models. The branching fraction 
i5(B — > xp{2S)X) is then derived from the fitted normalization. 

In addition to the fit uncertainty, including statistical and systematic uncertainties on the single 
measurements (±3.8%), the following sources of uncertainty are considered (with the corre- 
sponding relative AB uncertainty given in parentheses): 

• PDG branching fractions. The uncertainties quoted by the PDG for S(B ^ ]/tp X), 
B{xp{2S) ji^ji^), and B{]/tp H^f^^) are summed in quadrature (13.5%). 

• Model assumptions. The EvtGen prediction is used for an alternative fit and the 
difference with respect to the nominal N value is taken as a systematic uncertainty 
(1.0%). 

• FONLL uncertainties. All uncertainties on the underlying bb cross section, discussed 



in Section 8.1 are assumed to be fully correlated between B — > ]/xp X and B — > 
xp{2S)X transitions. Residual uncertainties affecting the cross-section ratio predic- 
tion are used to perform alternative fits and the differences with respect to the nom- 
inal value of N are taken as systematic uncertainties (1.5%). 

• B — > charmonium spectrum theoretical uncertainties. The only source of theoretical un- 
certainty which is not correlated is the one on the B — ?► charmonium spectrum. We 
use the high-Q and low-Q method to estimate this uncertainty as detailed in Sec- 
tioning In order to obtain an upper limit on the uncertainty of this ratio, the high-Q 
sample of B — ?► J/t/^X is compared to the low-Q B — ?► xp[2S)X sample and vice versa. 
The average difference in N with respect to the nominal EvtGen prediction is taken 
as a systematic uncertainty (3.8%). 

The measured value is: 

B{B tp{2S)X) = (3.08 ± 0.12(stat.+syst.) ± 0.13(theor.) ± 0.42(Spdg)) • 10^1 (12) 

where the last uncertainties are from the world-average branching fractions and the theoret- 
ical variations, respectively. The result is in agreement with the current world-average value 
from LEP and Tevatron measurements (6pdg(B ^ i/'(2S)X) = (4.8 ± 2.4) ■ 10"^ BSl), while 
improving the relative uncertainty by a factor of three. 



9 Summary 

A measurement of the ]/ip and ip{2S) production cross sections in pp collisions at -^s = 7TeV 
with the CMS experiment at the LHC has been presented. The data sample corresponds to 
an integrated luminosity of 36.7 ± 1.5 pb^^. The two cross sections and their ratio have been 
measured as a function of the meson transverse momentum (up to 70 GeV/c for the ]/ip and to 
30 GeV/c for the ip{2S)) in several rapidity ranges. Cross sections for prompt and non-prompt 
production have been determined from the measured values of the ]/tp and tp[2S) non-prompt 
fractions. 

The prompt cross section results are evaluated assuming isotropic decays in the production, as 
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well as four other polarization scenarios. In addition, cross sections restricted to the acceptance 
of the CMS detector are given, which are not affected by the polarization of the charmonium 

states. 

Cross sections for prompt and non-prompt production have been compared with NRQCD and 
FONLL predictions, respectively. Agreement is found in the prompt case: this is particularly re- 
markable in the ^{2S) case, where theoretical uncertainties are reduced because of the absence 
of feed-down from heavier charmonium states. In the non-prompt case, general agreement in 
shape is found for i/'(2S) in the entire pj range considered (up to 30 GeV/c), but an overall scale 
discrepancy is observed, possibly because of the assumption on the inclusive branching frac- 
tion B{B jp{2S)X). For J/1/7 there is similarly general agreement over the above range, while 
the predictions overestimate the measured differential cross-sections for 30 < px < 70 GeV/ c. 

For plausible hypotheses on the polarizations of the two charmonium states the ratio of their 
differential cross sections is obtained. In this ratio systematic errors largely cancel. The inclu- 
sive branching fraction B{B ^ xp(2S)X) is extracted from the ratio of the non-prompt cross 
sections to be: 

B{B ^p{2S)X) = (3.08 ± 0.12 (stat.+syst.) ± 0.13 (theor.) ± 0.42 (.Bpdg)) x 10^1 

improving the relative uncertainty on the previous world average by a factor of three. 
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